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■ Abstract. Buda-Lund hydro model fits are compared to BRAHMS, PHENIX, 
PHOBOS and STAR data on identified particle spectra, two-particle Bose-Einstein 
or HBT correlations, charged particle pseudorapidity distributions and pseudorapidity 

£SJ . as well as p t dependent elliptic flow in ^s NN = 130 and 200 GeV Au+Au collisions at 

\ RHIC. Preliminary results indicate that 7/8-th of the particle emitting volume is rather 

■ cold, with surface temperature of 105 MeV, but the temperature has a distribution 
and the most central 1/8-th of the volume is superheated to T{x) > T c — 172 ±3 MeV. 
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1. Introduction 

> ■ 

The Buda-Lund hydro model Ij is successful in describing the BRAHMS, PHENIX, 
PHOBOS and STAR data on identified single particle spectra and the transverse mass 
dependent Bose-Einstein or HBT radii as well as the pseudorapidity distribution of 
. £h ! charged particles in Au + Au collisions both at ^s NN = 130 GeV j2] and at a/s nn = 200 
GeV jS] • Recently, Fodor and Katz calculated the phase diagram of lattice QCD at finite 
net baryon density jlj. Their results, obtained with light quark masses four times heavier 
than the physical value, indicated that in the < //# < 300 MeV region the transition 
from confined to deconfined matter is not a first or second order phase-transition, but 
a cross-over with a nearly constant critical temperature, T c = 172 ± 3 MeV. The result 
of the Buda-Lund fits to Au+Au data, both at ^/s NN = 130 and 200 GeV, indicate 
the existence of a very hot region. The temperature distribution T(x) of this region 
is characterized with a central temperature To, found to be greater than the critical 
value calculated from lattice QCD: T > T c [3]. The Buda-Lund fits thus indicate quark 
deconfinement in Au + Au collisions at RHIC. The observation of a superheated center 
in Au+Au collisions at RHIC is confirmed by the analysis of p t and r\ dependence of the 
elliptic flow 0, measured by the PHENIX [Bj and PHOBOS collaborations [01 HDJ. A 
similar analysis of Pb+Pb collisions at CERN SPS energies yields central temperatures 
lower than the critical value, T < T c [HUT]. 

Here we summarize the Buda-Lund fit results as detailed in refs. El El IHj • 
See these papers for definitions and discussion of the results as well as more detailed 
references. 
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2. Buda-Lund fit results to central Au+Au data at ^s NN = 130 and 200 GeV 
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BudaLund v1 .5 fits to 200 AGeV Au+Au 
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Figure 1. The upper four panels show a simultaneous Buda-Lund fit to 0- 
5(6) % central Au+Au data at = 130 GeV, refs. QU d EE! El US] ' The 

lower four panels show similar fits to 0-30 % central Au+Au data at -\/s NN = 200 
GeV, refs. ^OEIEII- The fit parameters are summarized in Table 1. 

Analyzing the fit parameters (Table 1) we find that T > T c by more than 5cr in 
case of the — 5(6)% most central Au+Au data at ^/snn = 130 GeV. We interpret this 
as an indication of quark deconfinement. In case of the less central (0 — 30%) Au+Au 
data at ^snn = 200 GeV, with an improved analysis we find To > T c by 2a, not a 
significant difference. We interpret this as a possible hint for quark deconfinement. In 
both cases, the flow profile within errors coincides with the three-dimensional Hubble 
flow, u M = x^/t. For details, see refs. [HUH]- For similar results, see refs. (213 EH E2] • 
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3. Buda-Lund results for the elliptic flow in Au+Au data at RHIC 




Figure 2. In the top panels, the ellipsoidally symmetric Buda-Lund hydro calculations 
are compared to PHENIX identified v 2 {p t ) data at 200 GeV and PHOBOS u 2 (ry) 
data at 130 and 200 GeV 01113]. On the top right panel, we added a constant non- 
flow parameter of 0.02 to the calculated values of v 2 - In the lower panels, the flow 
coefficients, v n -s are shown for n — I ... 6 as a function of transverse momentum at 
midrapidity, as well as a function of pseudorapidity, with integrated pt ■ 

The model is evaluated in all panels with the following parameters: To = 210 MeV, 
T s = 105 MeV, X f = 8.6 fm, Y f = 10.5 fm, Z f = 17.5 fm, X f = 0.57, Y f = 0.45, 
Z f = 2.4, r = 7 fm/c, ■& = 0.09, ^ = 70 MeV, /i , K = 210 MeV and /i , P = 315 MeV. 
T is the temperature of the center, while T s is that of the surface, (Xf, Yf, Zf) are 
the principal axes and (Xf, Yf, Zf) the principal expansion rates of the ellipsoid at the 
freeze-out proper-time To. For more details, see pages 3-6 of ref. 

Table 1, Figures 1 and 2 indicate that the Buda-Lund hydro model works well 
at both RHIC energies: it gives a good quality description of the transverse mass 
dependence of the HBT radii as well as identified particle spectra and elliptic flow, 
see refs. [U EHl EH ES! for further details. 

The temperature was estimated to be above T c in the most central l/8th of the 
expanding ellipsoid [3], similarly to a fireball that is heated from inside. We interpret 
this result as a confirmation of the quark deconfinement and the cross-over like transition 
found in Buda-Lund fits to central Au+Au collisions at ^s NN = 130 and 200 GeV. 

This work was supported by following grants: OTKA T034269, T038406, the 
OTKA-MTA-NSF grant INT0089462, the NATO PST.CLG.980086 grant and the 
exchange program of the Hungarian and Polish Academy of Sciences. 
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Buda-Lund parameter 


Au+Au 200 GeV 


Au+Au 130 GeV 


71 [MpVI 

1 [1V1L V J 

T e [MeV] 
fin [MeV] 


117 ± 12 
61 ± 52 


91 A _|_ 7 

102 ± 11 
77 ± 38 


- n 'G L mi J 
ffml 

1^5 L J 

K) 


1 Q c + 17 

Id.tJ _I_ _L . 1 

12.4 ± 1.6 
1.6 ± 0.2 


8.6 ± 0.4 
1.0 ±0.1 


r [fm/c] 
At [fm/c] 


5.8 ± 0.3 
0.9 ± 1.2 
3.1 ± 0.1 


6.0 ± 0.2 
0.3 ± 1.2 
2.4 ± 0.1 


X 2 /NDF 


114 / 208 


158.2 / 180 



Table 1. Buda-Lund hydro model vl.5 source parameters, corresponding to Figs. 
1 and 2. The errors on the parameters are preliminary, as point-to-point and 
normalization errors are added in quadrature when evaluating \ 2 ■ 
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